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The characteristic formalism of numerical relativity is based on a system of coordinates aligned 
with outgoing null cones. While these coordinates were designed for studying gravitational waves, 
they can also be easily adapted to model cosmological past null cones (PNCs). Similar to obser- 
vational coordinates in the observational approach to cosmology, this then provides a model that 
only makes use of information causally connected to an observer. However, the diameter distance, 
which is used as a radial coordinate, limits the model's cosmological application to the region prior 
' to the PNC refocussing. This is because after refocussing, the diameter distance ceases to be a 

unique measure of distance. This paper addresses the problem by introducing a metric based on 
^ . the Bondi-Sachs metric where the radial coordinate is replaced by an afHne parameter. A model is 

Q ' derived from this metric and it is then shown how an existing numerical scheme can be adapted for 

\ simulation of cosmological PNC behaviour. Numerical calculations on this model are found to have 

. the same stability and convergence properties as the standard characteristic formalism. 

vn 

CN ' PACS numbers: 04.25.D-, 98.80.Jk 
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The observational approach to cosmology is the endeavour to reconstruct the geometry of the Uni- 
verse using observations that do not require the prior assumption of a cosmological model. Using this 
approach, it was shown in [l| that given ideal cosmological observations, the only essential assump- 
tion necessary to determine the geometry of the Universe is a theory of gravity. Assuming General 
I Relativity, then, the full set of Einstein field equations (EFEs) can be used to reconstruct the geom- 

. etry of the Universe using direct observations on the past null cone (PNC) as boundary conditions. 

\^ ' Observationally and theoretically this is a very ambitious task and therefore, current developments 

have been restricted to spherically symmetric dust models while only relaxing the usual assumption of 
homogeneity in the radial direction. These restricted models are important for the development of the- 
oretical foundations and also useful as verification models since they avoid the circularity of verifying 
what has already been assumed. For instance, doing investigations such as quantifying homogeneity 
on different scales, testing the verifiability of cosmology validating the Copernican principle 0] 
and determining the metric of the Universe Q require more general models than the conventional 
' Robertson- Walker geometries. 

5h ' 

' The Observational Cosmology (OC) programme of Ellis and others played an instrumental part 

in the development of the observational approach. Of particular interest is their use of observer 
coordinates which is based on the optical coordinates of Temple ^] . These coordinates are based on 
the natural propagation of electromagnetic radiation and follow the causal structure of an observer 
PNC. This approach then consists of two problems: firstly, astronomical observations are used to set 
up the metric on the local PNC and secondly these form the final values of a characteristic final value 
problem which determines the historical evolution of the region causally connected to the PNC (i.e. the 
interior of the PNC). The PNC and its interior, which is the causally connected region of a cosmological 
observer, is of fundamental importance since it defines the limits on which cosmological models can be 
validated from direct observations. Apart from the initial Ideal Observational Cosmology report 
which focussed on the fundamental principles, the focus of the OC programme has mostly been on exact 
solutions and perturbations using tetrad formalisms in spherical symmetry. Some recent developments 
are presented in @-(7'] and a summary in ^ . Numerical algorithms have been proposed in [3] using a 
metric based approach but these have not been implemented. A separate line of development based on 
the Lemaitre-Tolman-Bondi (LTB) model, using comoving coordinates, using the models introduced 
in , was followed in [ll|, [l^ to reconstruct the observer metric numerically. Some important 
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aspects addressed in the LTB approach are the treatment of a reconverging PNC through coordinate 
singularities and the sensitivity of the numerics when using reahstic data as boundary conditions. 

The characteristic formahsm of numerical relativity [isl - fisj usually associated with studying gravi- 
tational waves, was put forward in [l^[l3l as a method for modelling the observable universe. With the 
motivation that there exists a well-established base of null-cone numerical schemes, it was investigated 
to what extent these can be employed to simulate the evolution into the past of a cosmological PNC, 
instead of modelling the evolution into the future of a future null cone. It was found that with minor 
cosmological considerations, this method was very well suited for modelling the historical evolution 
of the observable universe. That is apart from a known limitation of the characteristic formalism; i.e. 
in modelling the reconvergence of the PNC of an expanding universe, the diameter distance, as radial 
coordinate, becomes multi-valued and limits the feasible radial extent of this formalism to the region 
prior to the observer apparent horizon (AH). Since the properties of the observable universe at the 
AH are directly related to the cosmological constant (A) d, [l^ , extending the model beyond the 
AH is a compulsory requirement for studying contemporary problems in cosmology. 

This paper addresses this problem by introducing a metric based on the Bondi-Sachs metric where 
the radial coordinate is replaced with an affine parameter. We derive the model with a cosmological 
constant A incorporated into the Einstein field equations, and regard A as a parameter of the theory of 
gravity rather than as a matter source term. Similar to the conventional characteristic formalism, this 
model consists of a system of differential equations for numerically evolving the EFEs as a characteristic 
initial value problem (CIVP). A numerical code implemented for the method has been found to be 
second order convergent. This code enables simulations of different models given identical data on the 
initial null cone and provides a method to investigate their physical consistency within the causally 
connected region of our current PNC. These developments closely follow existing 3D schemes developed 
for gravitational wave simulations, which should make it natural to extend the affine CIVP beyond 
spherical symmetric simulations. 

In section |TT1 an affine parameter is introduced in the Bondi-Sachs metric to derive the affine CIVP 
model for cosmology. Details of the numerical implementation are given in section IIIII and coordinate 
transformations, which relate comoving coordinates to the affine CIVP coordinates are discussed in 
section IIVI Results of simulations with the new coordinates are presented in section |Vl Section IVTl 
concludes the paper. 



II. CHARACTERISTIC FORMALISM ON NULL GEODESICS 

A. Characteristic formalism with affine radial coordinate 

The conventional characteristic formalism in numerical relativity uses a frame of reference based 
on outgoing null cones that evolve from values on an initial null cone. The idea is conceptualised in 
figure [T] G is a timelike geodesic, and u is the proper time on G. (In practice, in numerical relativity 
simulations the inner boundary of the null cone is usually a timelike worldtube rather than a geodesic, 
but the problem can be formulted with a geodesic). Null geodesies emanating from G have constant 
{u,9,ip), and near G the angular coordinates 9 and tp have the same meaning as in spherical polar 
coordinates. The coordinate r is the diameter, or area, distance defined by the condition that the 
surface area of a shell of constant r is Airr'^. 

The geometry of the characteristic formalism is described by the Bondi-Sachs metric, which in 
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spherical symmetry is formulated as ^ 

ds^ = ~e^f' (^1 + du" - 2e^l^dudr + r^d0^ + sin^ 0dip^}. (!) 

This can be recognised as a generalization of the well-known Eddington-Finkelstein form of the exterior 
Schwarzschild metric, obtained by setting /3 — 0,W — —2M where M is the mass of the source. The 
coordinate system is defined such that /3 and W vanish at the vertex of each null cone, i.e. at r = 
equation ([1]) reduces to a Minkowskian metric. 

Starting with the Bondi-Sachs metric ([1]), an affinely paramterized geodesic in these coordinates is 
determined through 

d^r 1 f dr\^ d'^r f dr\^ 

Setting A = r at the origin provides the initial conditions r(0) — A(0) — and dr/dX\x=o — 1 then 
solving gives 

Using tensor transformation laws and substituting ([3]) for all dr/dX terms, a new metric with the 
radial coordinate A is introduced as 

ds^ = -(^1 + y ^ du^ - 2dudX + f^{de^ + sin^ 0dip'^} (4) 

with: W = W{u,X) and r = f{u,X). 

Substituting (HJ into the EFEs, using the form Rat = n{Tab — \Tgab) + ^gab, with the stress-tensor 
for a dust-like fluid {Tab = pVaVb and T = —p) gives ^ 

-]j^i^rp{vif (5) 

r „A = ^ |w^,Ar A + ff^xx + Wr^xx - 2r „r a - 1 + (r a)^ + ^^pf^ + Af^j / f (6) 

W f 1 \ 

W^xx = — f,xx + 4r «A + 2k [vqVip - -p j f - 2Af (7) 

with: f(0) = W^(0) = M?,a(0) = r„(0) = and f a(0) = 1. 
Further, substituting the dust stress-tensor and (jj]) into the continuity equation, T°'\ — 0, the 



^ The notation used here is based on that of and substituting W = V — r will give the original notation of Bondi 

and Sachs in [H^. 
^ Using geometric units G = 1, c = 1 and k = Stt 
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energy-momentum equations follow 




with: Vw = 1 -\ — —■ 

r 

Making use of the normalisation condition, g°'^VaVb = —1, vq can be written in terms of vi as 

VO ^ ^Vn^Vl + ^V^\ (10) 

Having the values on the initial null cone for p and vi, equations ([5|) to ()10p form a hierarchical 
system that can be solved in the order (O, ([B]), ([7]) and (ITU)) , then solving equations ([S]) and © evolves 
the system to the next null cone where the process can be repeated until the domain of calculation has 
been covered. These equations are all interdependent and require an iterative scheme for a numerical 
solution. 



B. Cosmological considerations 

The model introduced in section fll Al ls essentially a null cone formalism which makes provision for 
a null cone that can reconverge at some distance from the cone vertex. Taking G to be the worldline of 
an observer located at the cone vertex and integrating into the past, these coordinates can be naturally 
aligned with a cosmological PNC. Further, it being spherically symmetric and radially inhomogeneous, 
classifies it as a LTB model in null coordinates. Using this model, the simulations done in and 
(l6j can be extended beyond the apparent horizon (AH) . Figure [1] illustrates the differences between 
conventional characteristic coordinates and the affine coordinates. Besides being not comoving, the 
affine coordinates closely resemble the observer coordinates as described in 




FIG. 1: On the left hand side, the null coordinates of the characteristic formalism and on the right, the afRne 
null cone coordinates for cosmology. 

As an observational cosmology problem, where the geometry is determined from direct observations, 
there are two subproblems to be solved: 



5 



i. Reconstruct the geometry of the local PNC from directly observable quantities, such as redshift- 
distance and galaxy number counts. 

ii. Evolve the model as a reverse CIVP to determine the interior of the PNC using the values 
determined in (i) as initial conditions. 

Although the accumulation of cosmological data in recent years has been astounding, these are 
not yet sufficiently complete for a practical implementation of the first problem, which will have to 
incorporate aspects such as data reduction and the sensitivity of the model to observational errors. 
These aspects will not be considered in this paper and only a conceptual description of the relations 
between observable quantities and the required input for the model will be described in the next 
section. The second subproblem will be treated in detail in the remainder of the paper through the 
development of a numerical code for the affine CIVP model. By itself this provides a mechanism for 
testing hypothesised models by evaluating the behaviour of their PNCs under given initial values. 
The PNC behaviours of different models given similar initial values are particularly interesting cases 
to investigate. 



C. Reconstructing the metric 



As initial data for the model, wi(A) and p{X) have to be determined from observations on the 
PNC. These require measures of the radial distribution of expansion and density. As a measure of 
expansion, redshift in terms of the luminosity distance (di) can be determined from redshift- magnitude 
observations (e.g. type la supernovae observations). The reciprocity theorem can then be used to 
convert di-redshift to diameter distance-redshift relations (see [22j) 

z{dL) and di, = {1 + z)'^ro ^ z = z{ro). (11) 

Here the zero subscript refers to a value on the current PNC. The redshift is directly related to the 
time component of the contravariant velocity: 

dU r, , , 

l + z= — =v° (12) 
dr 

where r is the proper time along general galactic world lines (see [l|). This can be used to determine 
the covariant velocity wi: 

v^{ro) = -v'' = -{l + z). (13) 

As a measure of density, observed galaxy number counts (n) in terms of z can be used to determine 
p{ro). Galaxy redshift surveys can provide data for this, although current surveys are not yet suffi- 
ciently complete at high redshifts. Since z(ro) is known, n(z) can be rewritten as n{ro) — r7,(z(ro)). 
It was shown in ^6|] that the proper number count can be written as 

In terms of the diameter distance, the proper number count then becomes 

Niro) - j^y (15) 
The proper density is then related to the proper number count 



p(ro) = f{N). 



(16) 
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The details of this relation will not be considered at this stage but in principle it must take into 
account aspects such as dark matter and source evolution, preferably with factors independent of an 
already assumed cosmological model. 

From ([Tn]) and fTC]) the required values of p and vi are obtainable in terms of the diameter distance. 
Since these are required in terms of A, interchanging the independent variable in ([5]) provides a method 
to relate A in terms of r. Starting with ([S]), rewritten as a system of ordinary differential equations 
(ODEs), introducing S, 

r A - S (17) 
S,\^ f^\\ = -]^nrp{viY (18) 
with: f (0) = and 5(0) = 1, 

then interchanging the roles of A and r, simplifying, rearranging and introducing U provides ODE 
expressions for \ f.f 

A,r = U (19) 
U^f - A,,^ = (}^nfp{vif^ il^ (20) 
with: A(0) = and U{Q) = 1. 



Although these equations will not be solved at this stage, there are some complications that will 
require special consideration around the AH. These arise from the fact that the diameter distance in an 
expanding universe is not necessarily monotonically increasing and therefore, not a unique independent 
variable. In these cases, f reaches a maximum and then decreases. At the maximum diameter distance, 
the AH, f.A = 0, and A.f is singular. In general, this can be overcome by separating the solution into 
regions prior, around and succeeding the AH where the region around the AH is solved using series 
expansions. Such a method was previously implemented in [llj to handle singularities around the AH 
for LTB models. 

It should be noted that Is^, followed by [6*1 and derived a relationship for the cosmological 
constant A (which we arc regarding as a parameter of the theory of gravity) involving the maximum 
value of f and matter data within the PNC. In order to test the theory of gravity more generally, 
a measurement involving a rate of change over time is required. Redshift drift z^u{z) 0], which is 
one of the design objectives of the CODEX spectrograph planned for the European Extremely Large 
Telescope (E-ELT) [23], can be used for this purpose. From (|13l) we have 

wi(ro) = -(1 + z)^ wi,«(?'o) -z^u{z), (21) 

and then an evaluation of the difference between the left and right hand sides of equation (|5]) provides 
a test of the theory being used. 

In the LTB models investigated in and [TS], data reduction methods were introduced which 
provide valuable insight into the methodology of converting realistic data into useful initial data. In 
combination with the methods introduced in this section, all options will have to be considered when 
data of sufficient completeness becomes available for observational models. 



III. NUMERICAL IMPLEMENTATION 



The numerical scheme described here is in structure similar to the finite difference scheme described 
in [l3|, which in turn was based on the general 3D code developed in in [l^ and [13]. Since the 
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original scheme is based on first order ODEs, the second order hypersurface equations are rewritten 
by introducing R{u, A) — f.„, S(u, A) and T{u, A), which then gives: 

r A = S (22) 
S,\ = r^AA = -^Kfp(i;i)^ (23) 

= f^„A = ^ l^fS + fS,x + WS^x -2RS-1 + S^ + ^Kpr^ + Af^j /f (24) 

Wx = f (25) 

T A = W^.AA = —5 a + 4 i?,A + 2« ivovip ~^P)f- 2Af (26) 
with: f(0) = M^(0) = i?(0) = f (0) = and S{0) = 1. 

These equations can be solved using standard ODE methods for systems of equations where (|22p and 
psp are solved as one system and ([24|) . (f25|) and (|26|) are solved as a separate system. In both cases 
a number of iterations are required to obtain convergence. 

The discretisation strategy followed is based on a rectangular grid similar to the one used in [l7j but 
using A as the radial coordinate. Solving the hypersurface equations is done with a central difference 
method on half steps between the r-grid points, using 

9j^9)-l + ^i9]rj+9"rj-l) (27) 

with i being the time step and j the radial step. Here j is calculated by substituting known values 
into equations (|22j|26l) . In order to solve the evolution equations, ([5]) and their general form is 
notated as 

ViM = F^i and = Fp (28) 
and as explicit finite differences on a time half step they are written as 



AuF:\^;'' and p]+' - p] + AuF;;''\ (29) 



Here, n is a time iterator that will approach i. In these equations, the numerical values at the point 
are used to evaluate the matter terms and hypersurface derivatives. Radial matter derivatives 
are calculated making use of standard central difference formulae (see for instance [1^ p. 160-161). 

After setting up a suitable grid, the numerical algorithm can be summarised in the following steps: 

i. Set the p and vi initial values on to the initial grid points. These values will, in principle, be 
obtained from observations. 

ii. Calculate r, r^u, W and va from p and vi on the initial null cone using (|22j|26p . 

iii. Calculate Fp, with n = 1 for the initial step, from the values of wi, p , r, r „, W and vq using 



1+1/2 ^ 

mation that will approach the actual values with subsequent iterations using (j29[) . 



IV. Set ^ Fp, again n = I for the initial step, and calculate vi and p as an initial approxi- 



V. Use the new values of vi and p to calculate r, r_„, W and vq and their radial derivatives similar 
to (ii). 

vi. Calculate ^"+^/^ = 1/2 (F" + F^'^^) from values in (v) and again vi and p for F^^^ and the 
values in (iii) for F^. 
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vii. Test the calculations in (vi) for accuracy and convergence. If they are sufficiently accurate, move 
to the next time step, otherwise repeat steps v. and vi. with the new values of vi and p. 



As with the standard characteristic model, described in [l6j_ 17 1, calculations in the regions around 
A = require special consideration. The mechanisms used in 17[ for these regions was directly adapted 
for the affine CIVP model. In the A « region, it is is evident from the occurrence of A denominators 
that equations ([S]) and © will not be well-behaved. Consequently, the A w region is calculated by 
making use of second order series expansions. The region where the series solution meets with the 
CIVP solution, also requires special treatment to avoid artificial instabilities. This has been done by 
smoothing out the merger region with a weighted average between the two solutions. 



IV. TRANSFORMATIONS: COMOVING COORDINATES 



A. General transformation 



As a measure of the degree of convergence and accuracy, numerical calculations have to be com- 
pared with results from known solutions, which require coordinate transformations from conventional 
comoving cosmological coordinates to affine null coordinates. These transformations follow directly 
from calculating the geodesic paths in space and time. This has the general form 

dx'' dx 



The parabolic LTB model, which is spherically symmetric and radially inhomogeneous will be used 
as the general case for transformations. Its geometry in comoving synchronous coordinates is described 
by the metric ^ 

ds^ = -df + [R,r{t, r)fdr^ + [R{t, r)]^{de^ + sin^ Od^^}. (31) 

Here, t is the cosmic (proper) time, r is a comoving radial coordinate with 6 and ip the inclination 
and azimuth angles. R{t, r) is the areal radius and AttR^ defines the proper surface area of a sphere 
with coordinate radius r at a constant time slice (26j . 

Using (pO)) . the geodesic equations for the LTB model becomes 



dH fdr 

dx^ + " Ua 

d^r _ „i dt dr / dr 



When scaled to some maximum time, to, as the current age of a universe, the conditions at A = are 
t ^ Iq and r = with the initial directions constrained to be null by dt/dX — 1 and dr/dX = i?(to, 0). 
As comparative values on a null cone grid, the covariant velocity follows directly from vi = dt/dX, 
the diameter distance from f(A) = R(t,r{X)), while p is determined from the coordinate expression 
for the specific model using (t(A). r(A)) as the coordinates for t and r on a null cone. 



^ Note that r here is the comoving radial distance and not the diameter distance. 
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B. Models for code verification 



1. Einstein- de Sitter model 



As an illustration of comoving to afBne transformations, the Einstein-de Sitter model (EdS), scaled 
to to = 1 with G = 1, which has the metric 

ds^ = -dt^ + t^l^dr^ + t^l^r^ide"^ + sin^ Odip^}. (34) 

provides a system of equations from which t{X) and r(A) can be obtained 

dH 2 , ,^ f drV 

d'^r A . dt dr „ 

dx^^rTxTx='- ^''^ 

At A = 0: t = 1, r = 0, dt/dX = -1 and dr/dX = t''^/^. 

Solving these equations numerically provides a useful example of the motivation for working with 
an affinely parameterized radial coordinate as opposed to the diameter distance. This is illustrated in 
figure [5] where the diameter distance and affine parameter is plotted against the redshift. In terms of 
observations, with the emphasis on the AH at z = 1.25, the diameter distance reaches its maximum 
and then desreases while A keeps on increasing and provides a unique coordinate for higher redshifts. 



0.6 
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Redshift (z) 



FIG. 2: The radial distance coordinates do and A related to z, which is a common observational measure. 

Substituting t{X) into the density equation p = l/(67rt^) provides the density profile on a specific 
null cone. It is interesting to point out here that as t varies with A on each null cone, p also varies i.e. 
even though the Einstein-de Sitter model is homogeneous in conventional cosmological coordinates, 
null cones are not hypersurfaces of radial homogeneity and the model becomes inhomogeneous in null 
coordinates (see [13] )• 



2. ACDM 



The ACDM or concordance model is of interest for the current understanding of the Universe and is 
also useful to test the affine CIVP model with a non-zero cosmological constant. In terms of the LTB 
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metric dSl]) . using cosmological properties at the current epoch (to), the solution of the flat ACDM 
model is (see [H]): 



i?(t,r) 



with the age of the Universe: 



and the density distribution 



1/3 



sinh 



-1 \ 2/3 



to = -[Ho\/^a] sinh" 



^rnO 



1/2 



(37) 



(38) 



3Hq rtjnoRo 



(39) 



Here, flmo and JIa are the current density parameters for Baryonic matter and the cosmological 
constant respectively and Hq is the current Hubble constant. Values representative of the actual 
Universe are: ^Imo = 0.3, J^a = 0.7 and Hq — 72 Mpc s~^ km~^. For the purpose of code verification, 
however, using geometric units (G = 1) and rescaling time to have to = 1, a dimensionless value of 
Ho = 0.964 was used to make comparison with other models more convenient. 



3. Lemaitre- Tolman-Bondi 



LTB models, introduced in |29l - l31j . have been popular in recent years to demonstrate how inho- 
mogeneities can reproduce the effects of type la supernovae redshift dimming without the need of a 
cosmological constant (see (3^-l36j amongst others). It, therefore, has physical significance and also 
provides a general case for testing the effect of radial inhomogeneity on the null cone behaviour. The 
test cases investigated, here, are for simplicity restricted to the parabolic solution. Substituting (j3ip 
into the EFEs and solving gives (see [iall^): 



^-M{r){t-ts{r)f 



and p(t,r) = — -f— . (40) 



47r R^R 



M{r) is the active gravitational mass which is the mass contributing to the gravitational field and 
tsir) is defined as the bang time function, which is a surface defined by the local time at which 
R = 0. A simplified subset of models with mathematical convenient properties follow by setting 
Mir) = Mq = r'^ as a coordinate condition where Mo is a constant which for illustrative purposes is 
set to 2/9. Equations (^0]) then reduce to: 

Rit,r)^r{t-tB{r))^^' (41) 

with p: 

''^ " 27r(i - tB{r)){3t - 3tB{r) - 2rtBAr)) ' ^^^^ 

By selecting t_B(r) = 0, equation (|^^ becomes the EdS model, which can be verified against the 
results in section HVB 11 If tB{r) ^ and iB,r(f) = 0, the time of the initial singularity is adjusted 
and the age of the Universe changes. For non-constant functions, the initial singularity becomes a 
singular surface and the age of the Universe becomes subject to the position of an observer (i.e. the 
age of the Universe depends on r). Thus, a variety of models can be generated for testing the code 
on parabolic spatial sections. 
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A simple choice of bang function is implemented as a verification model: 

tBir)=br, (43) 
with b being a constant. This simplifies equations (HIT) and (H^ to: 

R{t,r) ^r{t -hrfl^ (44) 

By varying the value of h, different aspects of inhomogeneity can be studied. The value 6 = is 
exactly the EdS model and 5 > shifts the age of a universe to a younger age as r increases while 
6 < provides the opposite effect where a universe is shifted to an older age as r increases. The latter 
case is particularly interesting since it provides a mechanism to mimic a cosmological constant on low 
redshifts [s^l- This is, however, not necessarily the case for high redshifts and the case where h = —0.5 
is used as a model where the effect of inhomogeneity is clearly visible on higher redshifts. 
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V. RESULTS 

A. Verification results 

The results of numerical calculations of the models described in section IIVBI are presented in this 
section. It is shown how the diameter distance (f), the covariant velocity (wi) and the density p evolve 
on PNCs in the past from observations on the current PNC. Also of interest is the evolution of f against 
the redshift z, which provides an illustration of an observer's perception of distance from a particular 
PNC vertex. In the results the current PNC values are indicated with uq, the oldest PNC with Umax 
and an intermediate PNC is also shown between these values. Since the purpose of the calculations is 
to test the behaviour of the numerical model and not the physics of the specific cosmological models, 
time is scaled to uq 1 a-nd geometric coordinates are used to provide dimensionless results. In all 
cases the numerical results (points in the figures) closely follow the transformed results (lines in the 
figures). Table U gives a summary of the test cases. 



Model 


Z 3,t Tfjiax 


Zmax 


uo 


"^max 


Figures 


EdS 


1.25 


4.5 


1 


0.35 


Emi 


ACDM [VIA = 0.7) " 


1.61 


6 


1 


0.4 


ElE] 


LTB (fo = -0.5) 


1.02 


1.2 


1 


0.2 


011] 



"This result is in agreement with that of [ssj l 

TABLE I: Affine CIVP test cases. 

Figure [7] B is a particularly interesting illustration of an observer's perception of distance in an 
inhomogeneous universe. Here, the (f, z) behaviour on earlier PNCs includes loops, which will com- 
pletely obscure the perception of diameter distance and redshift as measures of distance. This type 
of behaviour provides insight into the physical nature of a model by investigating its past behaviour. 
Although the initial PNC of the specific LTB model is not significantly out of line compared to more 
accepted models, especially on low redshifts, its past behaviour is rather unusual. This does not make 
it unrealistic, but does say that an observer at u = Umax would have found the interpretation of 
cosmological data particularly difficult. 



1. Einstein-de Sitter 
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(A) (B) 




FIG. 3: Diameter distance against A (A) and against of z (B) on PNCs at different proper times (w) evolved 
from a local PNC up to z = 4.5. 



(A) (B) 




FIG. 4: Density distribution (A) and covariant velocity (B) on PNCs at different proper times (u) evolved 
from a local PNC up to z = 4.5. 



2. KGDM with Q.a = 0.7 
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FIG. 5: Diameter distance against A (A) and against of z (B) on PNCs at different proper times (w) evolved 
from a local PNC up to z = 5. 
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FIG. 6: Density distribution (A) and covariant velocity (B) on PNCs at different proper times (u) evolved 
from a local PNC up to z = 5. 
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from a local PNC up to z = 1.2. 
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FIG. 8: Density distribution (A) and covariant velocity (B) on PNCs at different proper times (u) evolved 
from a local PNC up to z = 1.2. 

B. Convergence and accuracy 

Figures [9] A, B and C show the error propagation between the numerical calculations and the 
transformed results for different grid resolutions. Although the error values are generally small, it 
is visible that close to A = 0, where series expansions are used, the local error values are higher. 
No specific techniques were employed to reduce these errors, since this provides an indication of the 
model's stability to local errors. The convergence of errors are general and on higher grid resolutions 
the difference in error values are less evident. 

Figure [9] D, displays the convergence behaviour of the calculations. Comparing the error size of 
the p calculations with the grid size, gives the order of convergence between 1.7 and 2.2 i.e. around 
second order convergence. Since no specific techniques were used to balance the error sizes in the 
series expansion calculations with the CIVP calculations, the models also shows good convergence 
sensitivity against local errors. This effect will be investigated in detail in the future to determine 
how the model will behave with initial data with realistic observational error margins. 
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FIG. 9: Radial error distribution on the oldest null cones for the Einstein-de Sitter {umax = 0.35) (A), ACDM 
{umax = 0.4)(B) and LTB {umax ~ 0.2)(C) models for different grid resolutions. In (D), the error against grid 
size is shown; the slope of the change in error represents the order of convergence and is indicated in brackets 
in the legend. 
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VI. CONCLUSION 

The ideal observational approach proposed in [l| can arguably be considered the most philosophically 
sound approach to cosmology since it only studies the part of the Universe that is causally connected 
to a cosmological observer. Practically, a very precise map of the contents of the universe is required 
to develop the full potential of this approach, which, even in the era of precision cosmology, is not 
yet achievable. The more conventional model based approach is therefore, and will probably for 
some time, still be the most favourable approach for assessing our understanding of the Universe. 
However, the model and observational approaches do not have to be advocated in opposition since 
restricted implementations of the observational approach can provide insight into the question: Can 
our conventional understanding of the Universe, which is based on a priori assumptions such as 
the Copernican principle, be derived from less restrictive models making only the most essential 
assumptions? 

The afRne CIVP model developed in this paper provides a numerical implementation of the obser- 
vational approach restricted to a dust filled isotropic universe. The model is particularly well suited 
for investigating the historic evolution of the observable universe as a reversed CIVP evolved from 
observations on a local PNC. The line element of the affine CIVP, is closely related to the Bondi-Sachs 
metric, which makes it a direct process to adapt existing numerical schemes to the new metric. As 
with the standard characteristic formalism, second order convergence and accuracy was achieved when 
verified against models with known solutions. Similar to the simulations in [T7| . this provides a useful 
mechanism to investigate the physical consistency of cosmological models given identical observations 
on the observer PNC. With the standard characteristic formalism, this could only be implemented in 
the region prior to the AH, with the affine CIVP the region is extended far enough to investigate the 
effect of the cosmological constant on the history of the observable universe. 



Acknowledgments 

This work was supported by the National Research Foundation, South Africa. 



[1] G. F. R. Ellis, S. D. Nel, R. Maartens, W. R. Stoeger, and A. P. Whitman. Ideal observational cosmology. 

Phys. Reports, 124:315-417, 1985. 
[2] J.-P. Uzan, C. Clarkson, and G. F. R. Ellis. Time Drift of Cosmological Redshifts as a Test of the 

Copernican Principle. Physical Review Letters, 100(19):191303 — h. May 2008. 
[3] C. Hellaby and A. H. A. Alfedeel. Solving the observer metric. Phys. Rev. D , 79(4):043501, February 

2009. 

[4] G. Temple. New Systems of Normal Co-ordinates for Relativistic Optics. Royal Society of London 

Proceedings Series A, 168:122-148, October 1938. 
[5] M. E. Araiijo and W. R. Stoeger. Exact spherically symmetric dust solution of the field equations in 

observational coordinates with cosmological data functions. Phys. Rev. D , 60(10):104020, November 

1999. 

[6] M. E. Araiijo and W. R. Stoeger. Obtaining the time evolution for spherically symmetric Lemaitre- 

Tolman-Bondi models given data on our past light cone. Phys. Rev. D , 80(12):123517, December 2009. 
[7] M. E. Araujo and W. R. Stoeger. Using time drift of cosmological redshifts to find the mass-energy 

density of the Universe. Phys. Rev. D , 82(12):123513— 1-, December 2010. 
[8] M. E. Araujo and W. R. Stoeger. Finding a spherically symmetric cosmology from observations in 

observational coordinates - advantages and challenges. JCAP, 7:29 — h, July 2011. 
[9] N. Mustapha, C. Hellaby, and G. F. R. Ellis. Large-scale inhomogeneity versus source evolution - Can 

we distinguish them observationally? Mon. Not. R. Astron. Soc, 292:817 — h, December 1997. 
[10] N. Mustapha, B. A. C. C. Bassett, C. Hellaby, and G. F. R. Ellis. The distortion of the area distance- 

redshift relation in inhomogeneous isotropic universes. Classical and Quantum Gravity, 15:2363-2379, 



18 



August 1998. 

[11] T. Hui-Ching Lu and C. Hellaby. Obtaining the spacetime metric from cosmological observations. Classical 
and Quantum Gravity, 24:4107-4131, August 2007. 

[12] M. L. McClure and C. Hellaby. Determining the metric of the Cosmos: Stability, accuracy, and consis- 
tency. Phys. Rev. D , 78(4):044005, August 2008. 

[13] N. T. Bishop, R. Gomez, L. Lehner, M. Maharaj, and J. Winicour. High-powered gravitational news. 
Phys. Rev. D , 56:6298-6309, November 1997. 

[14] J. Winicour. Characteristic Evolution and Matching. Living Reviews in Relativity, 12:3, April 2009. 

[15] C. Reisswig, N. T. Bishop, D. PoUney, and B. Szilagyi. Unambiguous determination of gravitational 
waveforms from binary black hole mergers. Phys. Rev. Lett., 103:221101-1-221101-4, 2009. 

[16] N.T. Bishop and P. Haines. Observational Cosmology and Numerical Relativity. Quaest. Math., 19:259- 
274, 1996. 

[17] P. J. van der Walt and N. T. Bishop. Observational cosmology using characteristic numerical relativity. 

Phys. Rev. D, 82(8):084001, Oct 2010. 
[18] C. Hellaby. The mass of the cosmos. Mon. Not. R. Astron. Soc., 370:239-244, July 2006. 
[19] H. Bondi. Gravitational waves in general relativity. Nature, 186:535-535, 1960. 

[20] H. Bondi, M. G. J. van der Burg, and A. W. K. Metzner. Gravitational Waves in General Relativity. VH. 
Waves from Axi-Symmetric Isolated Systems. Royal Society of London Proceedings Series A, 269:21-52, 
August 1962. 

[21] R. K. Sachs. Gravitational Waves in General Relativity. VHI. Waves in Asymptotically Flat Space- Time. 

Royal Society of London Proceedings Series A, 270:103-126, October 1962. 
[22] G.F.R Ellis. Relativistic Cosmology. In General Relativity and Gosmology, Proc. Int. School of Physics 

'Enrico Fermi' (Varenna), Course XLVII, Ed. R.K. Sachs, pages 104-179. Academic Press, 1971. 

Reprinted in Gen. Rel. Grav. 41, 581 (2009). 
[23] L. Pasquini, S. Cristiani, R. Garcia-Lopez, M. Haehnelt, and M. Mayor. CODEX: An Ultra-stable High 

Resolution Spectrograph for the E-ELT. The Messenger, 140:20-21, June 2010. 
[24] N. T. Bishop, R. Gomez, L. Lehner, M. Maharaj, and J. Winicour. Incorporation of matter into charac- 
teristic numerical relativity. Phys. Rev. D , 60(2):024005, July 1999. 
[25] R. L. Burden and J.D. Faires. Numerical Analysis. PWS Publishing Company Boston, fifth edition, 1993. 
[26] G. F. R. Ellis and H. van Elst. Cosmological Models (Cargese lectures 1998). In NATO ASIC Proc. 541: 

Theoretical and Observational Cosmology, pages 1-116, 1999. 
[27] M. B. Ribeiro. Observations in the Einstein-De Sitter cosmology: Dust statistics and limits of apparent 

homogeneity. Astrophys. J. , 441:477-487, March 1995. 
[28] J.E. Lidsey. ASTM108 Cosmology/MTH703U Advanced Cosmology. Lecturer's notes. 

http://www.maths.qmul.ac.Uk/~jel/ASTM108/#coursenotes (cited: January 2010). 
[29] G. Lemaitre. I'Univers en expansion. Annales de la Societe Scientifique de Bruxelles, 53:51 — h, 1933. For 

an English transaltion see: The Expanding Universe. Gen. Rel. Grav, 29(5):641, 1997. 
[30] R.C. Tolman. Effect of inhomogeneity on cosmological models. Proc. Nat. Acad., Sci. 20:169, 1934. 
[31] H. Bondi. Spherically symmetrical models in general relativity. Mon. Not. Roy. Astr. Soc, 107:410, 1947. 
[32] M.-N. Celerier. Do we really see a cosmological constant in the supernovae data? A & A, 353:63-71, 

January 2000. 

[33] J.-F. Pascual-Sanchez. Cosmic Acceleration:. Inhomogeneity Versus Vacuum Energy. Modern Physics 
Letters A, 14:1539-1544, 1999. 

[34] K. Enqvist. Lemaitre Tolman Bondi model and accelerating expansion. General Relativity and Gravita- 
tion, 40:451-466, February 2008. 

[35] J. Garcia-Bellido and T. Haugb0lle. Confronting Lemaitre Tolman Bondi models with observational 
cosmology. JCAP, 4:3, April 2008. 

[36] M.-N. Celerier. The Accelerated Expansion of the Universe Challenged by an Effect of the Inhomogeneities. 
A Review. arXiv:astro-ph/0702416, February 2007. 

[37] J. Plebahski and A. Krasihski. An Introduction to General Relativity and Cosmology. Cambridge Univer- 
sity Press, first edition, 2006. 

[38] M. E. Araujo and W. R. Stoeger. The angular-diameter distance maximum and its redshift as constraints 
on A / Friedmann-Lemaitre- Robertson- Walker models. Mon. Not. Roy. Astr. Soc, 394:438-442, March 
2009. 



